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ABSTRACT 

Aims. We present Subaru/Suprime-Cam deep V and / imaging of seven fields in the outer regions of M33. Our aim is to search for 
stellar structures corresponding to extended Hi clouds found in a recent 21 -cm survey of the galaxy. Three fields probe a large Hi 
complex to the southeastern (SE) side of the galaxy. An additional three fields cover the northwestern (NW) side of the galaxy along 
the Hi warp. A final target field was chosen further north, at a projected distance of approximately 25 kpc, to study part of the large 
stellar plume recently discovered around M33. 

Methods. We analyse the stellar population at R > 10 kpc by means of V, I colour magnitude diagrams reaching the red clump. 
We constrain the age and metallicity of the different stellar populations, search for density enhancements that correspond to the Hi 
features, and investigate the radial surface distribution of the stars. 

Results. We find evolved stellar populations in all fields out to 120' (~ 30 kpc), while a diffuse population of young stars (~ 200 
Myr) is detected out to a galactocentric radius of 15 kpc. The mean metallicity in the southern fields remains approximately constant 
at [M/H]= -0.7 beyond the edge of the optical disc, from 40' out to 80'. Along the northern fields probing the outer Hi disc, we also 
find a metallicity of [M/H]=-0.7 between 35' and 70' from the centre, which decreases to [M/H]= -1.0 at larger angular radii out to 
120'. In the northernmost field, outside the disc extent, the stellar population of the large stellar plume possibly related to a M33-M31 
interaction is on average more metal-poor ([M/H]=-1.3) and older (> 6 Gyr). 

Conclusions. An exponential disc with a large scale-length (~ 7 kpc) fits well the average distribution of stars detected in both the 
SE and NW regions from a galactocentric distance of 1 1 kpc out to 30 kpc. The stellar disc extends beyond the Hi disc. The stellar 
distribution at large radii is disturbed and, although there is no clear correlation between the stellar substructures and the location of 
the Hi clouds, this gives evidence of tidal interaction or accretion events. 

Key words. Galaxies: M33 - Galaxies: evolution - Galaxies: stellar content - Galaxies: halos 



1. Introduction 

Hierarchical models of structure formation predict that galax- 
ies form after a series of accretion events involving lower mass 
systems ( |Abadi et al.|2003||Governato et al.|20Q7| >. Within this 
scenario, a massive spiral galaxy such as the Milky Way should 
have accreted 100-200 satellites during its evolutionary history 
(Bullock & Johnston 2005). This process would leave its im- 
print in a diffuse stellar halo extending out to 10-100 times the 
optical radius of the galactic disc. The density, luminosity, and 
metallicity of the stellar halo can thus provide direct information 
about the history and evolution of a galaxy, and enable us to test 
hierarchical formation scenarios. 

The discovery of the Sagittarius stream in the Milky Way 
(MW) revealed for the first time the role of satellite accretion 
in the build-up of a galaxy mass (Ibat a et aI7|| 1994| >. Studies of 
the MW halo using data from the Sloan Digitised Sky Survey 
(SDSS; |York et aT7||2QOO| > revealed other tidal streams around 
our galaxy, proving that these features are relatively common 
and can survive for a few Gyrs ( Belokurov et al. 2006| Grillmair 



[20061 |Belokurov et al.|2007)|Grillmair|2 009 ). 



* Based on data collected at Subaru Telescope, which is operated by 
the National Astronomical Observatory of Japan. 



Similar features were discovered around M3 1, the other mas- 
sive spiral of the Local Group (LG), that implied this galaxy 
has experienced a more violent merger history than the MW 
( |Font et al.|2006[|Fardal et al.|200"8"] >. A giant stellar stream with 
a neutral hydrogen (Hi) counterpart was found in this galaxy 
(|Newton & Emerson||1977| |Ibata et aL]|2001| |Ferguson et al.| 
|2002[|Thilker et al.|2004| >. Much effort has been devoted to ob- 
taining deeper optical images of the halo of M3 1 with the Pan- 
Andromeda Archaeological Survey (PAndAS; |McConn achie 
|et al.|2009) , which revealed a large-scale metal-poor stellar dis- 
tribution (stretching out to 150 kpc), and uncovered additional 
extended substructures in the halo similar to the giant southern 
stream. 

While M3 1 clearly shows signs of disturbance in its spheroid 
and outer disc, M33, its brightest satellite and the third most 
massive member of the LG, has long been considered an exam- 
ple of a "pure disc" system evolving in relative isolation, because 
of the lack of a prominent bulge, a seemingly unperturbed stellar 
disc, and the absence of known nearby satellites. The origin of 
bulgeless galaxies such as M33 is difficult to explain within the 
hierarchical structure formation scenario ( |van den Bosch|200T] l, 
because mergers of galaxies are expected to form a bulge compo- 
nent. The inflow of gas at later times (z > 1) from the local envi- 
ronment is thought to play a major role in the growth of the discs 
of these systems ( |Dekel & Bknboim|2006} |Brooks et al.|2009) >. 
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Nonetheless, the presence of a diffuse faint stellar halo even in 
disc galaxies seems to be expected within the framework of cur- 
rent galaxy formation scenarios, predicting that a diffuse stellar 
halo rather than a centrally-concentrated bulge would result af- 
ter the late major-merger event experienced by these galaxies 
dBrook et al.|20T0) . 

Studies of the outskirts of M33 have started providing evi- 
dence of an extended stellar distribution beyond the optical disc, 
although it remains unclear whether this component corresponds 
to a stellar halo or an extended disc. |Brooks et al.| ( |2004[ ) derived 
a radial density stellar profile out to 1 degree from the centre of 
M33 finding a clear dropoff in the star counts as a function of 
radius (oc R l - 46 ), that they interpreted as a evidence of a stel- 
lar spheroid component. The metallicity peak of this sparse stel- 
lar population, [Fe/H] = -1.24 + 0.04, was lower than the abun- 
dance in the disc ([Fe/H] — 0.9), and comparable to the mean 



metallicity of a sample of halo globular clusters (Sarajedini et al. 



2000 ). In a spectroscopic survey of red giant branch (RGB) stars 
in M33, |McConnachie et al.| ( |2006[ ) found that the stellar ve- 
locity distributions could be described by three different com- 
ponents: the galaxy disc, a halo (with [Fe/H] = -1.5), and a third 
distribution possibly associated with a stellar stream. Finally, the 
PAndAS survey revealed an extended stellar distribution around 
M33 well beyond the main optical disc out to 3 degrees (or 
a projected radius of ~ 45 kpc) from its centre ( McConnachie 
|et al |2009||2010| l. This structure could either be the stellar de- 
bris of an accreted satellite, or the result of a tidal interaction 
with M3 1 . This latter interpretation would be supported by the 
similar spatial orientation of the stellar features and the M33 Hi 
warp. The origin of the Hi warp remains disputed and previ- 
ous studies have proposed that it is the result of an interaction 
with its massive neighbour ( Corbelli & Schneider|1997} [Bekki 
|2008|l. Simulations of the relative motion of the M31/M33 sys- 
tem (|Loeb et al.|2005)|Putman et al.|2009[|McConnachie et al.| 
|2009| l, converge on the possibility of a close encounter between 
the two galaxies occurring a few Gyr ago with a pericenter dis- 
tance around or greater than 40 kpc, which could strip material 
from the M33 outer disc without severely distorting or disrupting 
the bulk of the stellar disc. 

Additional information about the outer regions of M33 can 
be inferred from neutral hydrogen studies. A 21 -cm survey of the 
Hi distribution in a region of 3°x3° around M33 with the Arecibo 
telescope (ALFALFA; Giova nelli et al.|2003) led to the finding 
of a popul ation of Hi clouds in the environment of M33 (Grossi 
et al.|2008 1. These clouds, with Hi masses rangi ng between 10 



and few times 10 6 M (at a distance of 840 kpc; |Freedman et al. 
2001), were found within a projected distance of about 20 kpc, 
and they appear to be distributed (in projection) along the major 
axis of the Hi disc, towards the direction of M3 1 . These features 
could also be evidence of either a tidal interaction with M3 1 , or 
the accretion of a dwarf companion of M33. In either case one 
might expect to find features in the stellar distribution that corre- 
spond to these Hi structures. The search for optical counterparts 
to the most massive and extended clouds may provide further in- 
sight into their origin, and the understanding of the outskirts of 
M33. 

Taking advantage of the unique sensitivity and large field 
of view of the Subaru Suprime-Cam instrument we mapped 
through the V and / filters sixregions in the southeastern (SE) 
and northwestern (NW) outer disc of M33 corresponding to the 
main Hi structures found in Gros si et al.| ( |2008"] l. An additional 
field to the north of the galaxy was observed to study in more de- 
tail the nature of the stellar structure at very large radii found by 
|McConnachie et aT7| ( [2009| l. We used colour-magnitude diagrams 
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Fig. 1. Location of the seven target fields around M33 imaged 
with Subaru/Suprime-Cam. The M33 Hi disc is highlighted by 
two contours with a column density of N#/ = 5 and 50 x 10 19 
cm -2 , respectively. Contours are overlaid on the GALEX (far 
ultraviolet) image of M33 ( |Thilker et aLl|2005b| >. The main Hi 
features to the SE and NW of the galaxy are also shown. 



(CMDs) reaching down to the red clump (RC) to constrain the 
age and the metallicity of the stars at large galactocentric dis- 
tances, and we searched for stellar density enhancements related 
to the gaseous features. 

The paper is organised as follows: in Sect. 2, we describe the 
observations and data reduction process; in Sect. 3, we discuss 
the main properties of the CMDs of the target fields; in Sect. 4, 
we derive the distance to M33 from the tip of the RGB. We esti- 
mate the age and metallicity of the stellar populations in Sect. 5 
from the comparison to evolutionary stellar population models, 
we derive in Sect. 6 the star formation history of the northern- 
most field of our survey using the CMD fitting method, while in 
Sect. 7 we show the spatial distribution of the stellar structures 
in the target regions. Discussion of our results is given in Sect. 
8, and in Sect. 9 we summarize our conclusions. 



2. Observations and data reduction 



We obtained V and / images with Suprime-Cam (Miyazaki et al. 
|2002[ ) on the 8.2 meter Subaru telescope on the nights of August 
21 and September 21-22, 2009. Suprime-Cam consists of ten 
CCDs of 2048x4096 pixels with a pixel size of 0.2 arcsec/px 
and a total field of view of 34x27 arcmin. The regions observed 
during the run are shown in Fig. [T] The fields are labelled from 
one to seven and are related to various Hi features around the 
galaxy. Fields SE1-SE3 cover the large gaseous complex to the 
south-east of the galaxy (see Fig.l), whereas fields NW4-NW6 
extend from the northern edge of the disc along the Hi warp to 



the northernmost Hi clouds detected around the galaxy in Grossi 
|et al.| ( [2008] l. Field NW7 was chosen to analyse the stellar pop- 
ulation of the large stellar feature uncovered by the PAndAS 
Survey. Finally, a field at approximately 4 degrees to the east 
of the optical center of M33 was selected to inspect the contami- 
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Table 1. Central coordinates of the target fields observed with 
Subaru/Suprime-Cam, corresponding projected distance from 
M33 (for dM33 = 840 kpc), seeing of the combined V and / im- 
ages, and extinction values adopted for each field. 



Field 


RA 


DEC 


r 


Seeing 


A v 




(J2000) 


(J2000) 


kpc 


V 1 I 


mag 


SE1 


01:36:15.6 


29:40:30 


16.3 


0".9/0".8 


0.157 


SE2 


01:36:51.6 


30:06:00 


12.6 


l".l/0".8 


0.179 


SE3 


01:36:22.8 


30:31:30 


8.2 


0".7/0".6 


0.184 


NW4 


01:32:37.6 


31:29:24 


12.7 


0".6/0".6 


0.151 


NW5 


01:30:28.8 


31:34:12 


17.0 


0".7/0".6 


0.177 


NW6 


01:28:20.0 


31:39:00 


22.6 


0".7/0".6 


0.152 


NW7 


01:31:46.9 


32:20:00 


25.4 


0".6/0".6 


0.176 


BG 


01:50:50.9 


30:39:36 




0".6/0".6 


0.168 



nation from foreground and background sources. For each field, 
we obtained a set of 5x440s and 15x320s exposures through the 
V and / filters, with total exposure times of 2200s and 4800s, re- 
spectively. To fill the gaps between the CCD chips and improve 
the removal of the effects of cosmic rays and bad pixels, the indi- 
vidual exposures were dithered by 60 arcsec along a five position 
pattern. Table 1 lists the coordinates of the eight target fields and 
the corresponding projected distance from the centre of M33. 



2.1. Photometry 

The images were bias-subtracted, fiat-fielded, distortion- 
corrected, and combined using the standard procedures of the 
Suprime-Cam Deep Field Reduction (SDFRED) software ( |Yagi| 
et al.|2002 Ouchi et al.|2004 1. The photometry was performed 
using the IRAF version of the package DAOPHOT ( |Stetson|1987) . 
A preliminary selection of detections with signal-to-noise ratio 
greater than 4 was performed with the automatic star-finding 
algorithm DAOFIND and the aperture photometry was obtained 
with the task PHOT using an aperture radius of 9 pixels for the 
images taken during the August run (fields SE1 and SE2), which 
had the poorest seeing (see Table 1), and an aperture of 6 pixels 
for the others. An empirical point-spread function was built with 
the task PSF by selecting several isolated and bright stars in each 
field. Finally PSF-fitting photometry was extracted with the task 
ALLSTAR (Stetson 1994). To reduce spurious detections, the fi- 
nal catalogs for each field includes only objects with^ 2 < 3, and 
| sharp | < 0.3, where x 1 is a parameter that indicates the quality 
of the fit, and sharp measures the spatial extent of a detection, 
being defined as the difference between the square of the width 
of the measure object and the square of the width of the PSF. 

The photometry was calibrated using data for Landolt ( |1992| l 
standards observed during the entire run. Galactic extinction cor- 
rections were derived following Schl egel et al.| ( fT9"98 i. We ap- 
plied an average correction corresponding to the extinction at 
the central position of each field. The adopted values of Ay are 
displayed in Table 1. 



2.2. Completeness 

We performed artificial star tests to evaluate the completeness of 
our photometry. We added artificial stars to the original images 
with magnitudes in the ranges 19 mag < / < 26 mag, and 19 mag 
< V < 27 mag, with a binning step of 0.2 mag that was increased 
to 0.5 mag for magnitudes brighter than I = 23 mag. A total of 
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Fig. 2. CMDs of fields SE1 (top-left), SE2 (top-right), SE3a 
(bottom-right), and SE3b (bottom-right). The boxes in the top- 
left panel display the location of RGB and RC stars in the dia- 
grams. Ten Gyr isochrones from Mari go et al.| ( |2008| at a metal- 
licity of [M/H]= -0.4, -0.7, -1, -1.3 are overlaid on the top-right 
panel to give an indication of the expected metallicity range of 
the RGB population. Light blue and red tracks in field SE3a 
(bottom-left) correspond to a metallicity [M/H]= -0.7 and an age 
of 250 Myr, 1 Gyr, and 10 Gyr, while dark blue and purple tracks 
in field SE3b (bottom-right) indicate the same age range for a 
metallicity of [M/H]= -0.4. 



56000 stars were injected at random positions onto each image 
for both filters. Photometry was then performed in the same way 
as described in the previous section, and the same selection cri- 
teria on the DAOPHOT parameters were applied. The final de- 
tections were considered as stars if their position was within two 
pixels of the original list of injected sources, and their magni- 
tude within 0.3 mag of the assigned value. Our photometric data 
reach the 50% completeness level at V ~ 25.6 mag and / ~ 24.7 
mag in all fields. 



3. Colour-magnitude diagrams 

3.1. The southeastern fields in the region of the Hi complex 

Figure [2] shows the CMD of the three fields in the southeast- 
ern region of the galaxy where the main Hi complex is located 
(see Fig. 1). The CMDs of field SE1 and SE2, which cover the 
two main substructures of the southeastern gaseous complex, are 
displayed in the top-left and top-right panels, respectively. Field 
SE3 maps the region between the Hi feature and the optical disc 
of M33; the CMD of this area is shown in the bottom panels of 
the figure. To differentiate the stars at the edge of the optical disc 
from a more diffuse stellar component, we divided the field into 
two halves (SE3a, and SE3b; see Fig. 1) and display the corre- 
sponding diagrams in the bottom-left and bottom-right panel of 
the figure, respectively. 
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Fig. 3. CMDs of field NW4 (top-left), NW5 (top-right ), NW6 
(bottom-left), and NW7 (bottom-right). Isochrones from Marigo 
|et al.| ( |2Q08| > are overlaid on the diagrams. Blue and red tracks 
in the first three panels indicate a metallicity of [M/H]= -0.7 
and ages of 250 Myr, 1 Gyr, and 10 Gyr. The same metallicity 
and age range is shown in the diagrams to mark the change in 
the RGB shape from NW4 to NW6. Isochrones in field NW7 
(bottom-right) correspond to an age of 10 Gyr and a metallicity 
of [M/H]= -0.7, -1, -1.7 (fromright to left) and give an indication 
of the expected metallicity range in the RGB stellar population. 



The most prominent feature is the red giant branch (RGB), 
corresponding to stars with colours (V - 7)o between 0.7 and 
2 mag, and Iq < 24 (see top-left panel in Fig. [2j. The RGB is 
clearly detected in all fields, although in SE2 and SE3a this fea- 
ture is relatively less populated, especially in the upper end. The 
edge of the optical disc of the galaxy is shown in panel SE3b as 
indicated by the redder and more populated RGB. 

Because of the well-known age-metallicity degeneracy, dif- 
ferent combinations of age and metal abundances are consistent 
with the colour of this feature. This is shown in the top-right 
panel where the RGB of SE2 is compared to isochrones from 
Marigo et al. (2008) with a common age of 10 Gyr and metallic - 
ities of [M/H]= -0.4, -0.7, -1,-1.3 (from redder to blue colours), 
where [M/H]= log(Z/Z Q ) and Z Q = 0.019. On the other hand, 
the width of the RGB could also be due to an age spread as the 
two isochrones at 1 and 10 Gyr (with [M/H]= -0.7 in SE3a, and 
[M/H]=-0.4 in S3b) show. In particular for SE3b, the age range 
of the stellar populations may be several Gyrs. 

Figure[2]also indicates a population of stars in the red clump 
(RC) phase, i.e. intermediate-age (1-10 Gyr) stars experienc- 
ing a core-helium burning phase. The average magnitude of 
this feature is around Iq ~ 24.4 mag, and the colour range is 
0.7 < (V - I)o < 1.2 (as shown by the box in the top-left panel 
of Fig. [2}. In addition the RC mean colour and magnitude de- 
pend on the age and metallicity of the stars ( |Girardi & S alaris 
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Fig. 4. CMD of the control field used to inspect the foreground 
and background contamination. 



for stars older than 2 Gyr, is mainly sensitive to metallicity vari- 
ations (Girardi & Salaris 2001). We use these properties of the 
RC feature to constrain the age and metallicity of these stars in 
the different fields (see Sect. 5). 

The CMDs of SE3 show a small population of main- 
sequence blue stars with (V - I)o < that are both related to 
the edge of the optical disc (SE3b), and scattered throughout the 
field out to a maximum projected distance of 50' (SE3a), about 
ten times the disc scale (6'4 in the J band; Regan & Vogel 1994 1. 
Comparing the brightest magnitude of the blue stars (i.e. the 
main-sequence turn-off luminosity) to isochrones with [M/H]=- 
0.7 implies an age around 250 Myr. 



2001 ), although the magnitude of the RC is more sensitive to 



age rather than to metallicity ( jDurrell et al.|2010"| , and its colour, 



3.2. The northwestern fields 

3.2.1 . Stellar populations in the Hi warp 

Fields NW4-NW6 map the northwestern side of the galaxy along 
the Hi warp out to the northernmost Hi clouds found around 
M33. The CMDs of field NW4 (Fig. [3] top-left panel) show a 
stellar population with similar properties as the southern regions 
at a comparable distance from the centre of M33 (SE2). On the 
other hand, the RGB of NW5 and NW6 (Fig. [3] top-right and 
bottom-left panels) gradually becomes bluer, narrower, and less 
populated, implying a smaller age spread in the RGB, and/or a 
change in the average metal abundance. Isochrones at 1 and 10 
Gyr with [M/H]= -0.7 are overlaid as a comparison. A small 
population of blue main-sequence stars is also detected in NW4 
and NW5, out to a maximum angular distance of 60' . 



3.2.2. Stellar populations in the large stellar feature around 
M33:the NW7 field. 

Field NW7 is not related to any Hi feature since this region was 
not included in our 21-cm survey of M33, but was selected to 
investigate the stellar population of the puzzling large stellar 
feature around M33 found by the PAndAS survey. It is roughly 
one degree further north in declination than NW6. Despite be- 
ing at approximately the same angular distance from the centre 
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Fig. 6. CMDs of southern fields in radial bins of width 10'. The radial angular distance from the M33 centre is displayed at the 
top-left corner of each plot. The isochrone that best fits the RGB shape and colour is also shown. The corresponding metallicity and 
age are indicated in the top-left corner of each panel. 
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Fig. 5. Contamination-subtracted Hess diagrams of field SE1 
(top-left), SE2 (top-right), NW5 (bottom-left), and NW7 
(bottom-right). The resolution of the diagrams is 0.05 x 0.05 
mag, and the colour scale indicates the number of stars per pixel. 
Only positive residuals are shown in the diagrams. 



of M33, the difference between the CMD of these two fields is 
striking. NW7 shows a clearly detected RGB, and RC, in con- 
trast to NW6, where the RC is hardly visible. While the number 
of RGB stars decreases from NW4 to NW6, NW7 shows a pe- 
culiar increase in the stellar density, given that the correspond- 
ing projected radius is around 25 kpc (for a M33 distance of 840 
kpc). Moreover, the overall colour of the RGB and RC features is 
bluer than in the other fields as can be inferred from the compar- 



ison to the stellar isochrones. Population synthesis models with 
a common age of 10 Gyr and [M/H]= -0.7, -1, -1.7 are over- 
laid on the CMD in the bottom-right panel of Fig. [3] They show 
that the stellar metal abundances could range between [M/H]=- 
1.7 and [M/H]=-l, suggesting the presence of a population with 
different properties from the other fields. 



3.3. Foreground and background contamination 

We observed a field offset by 4 degrees from the centre of M33 to 
take into account the foreground and background contamination 
due to Milky Way stars and unresolved distant galaxies misclas- 
sified as point sources. The field was processed in the same way 
as the others and the same cuts were applied to the final list of 
detections. The CMD of this region is displayed in Fig. [4] and 
shows that contamination mainly occurs in the region between 
0.3 < (V - I)o < 0.8 and Iq > 24 mag, and at colours redder 
than (V - I)o ~ 1.5 mag. Foreground and background contam- 
inants may affect the lower part of the RC, while the region of 
the RGB and blue main sequence are not severely contaminated. 
At (V - I)o < 1 mag and Iq > 23 mag, the main sequence of 
stars in the Galactic halo is clearly visible, but mainly contami- 
nates the blue and bright region of the CMDs, where we do not 
find significant stellar populations related to M33. To assess the 
contribution from foreground and background contaminants, we 
built the Hess diagram of the control field, providing the den- 
sity of stars at different positions in the CMD, dividing the di- 
agram into bins of dimension 0.05 x0.05 mag in colour and / 
magnitude. We then subtracted it from the Hess diagrams of the 
other fields around M33, and, as an example, in Fig. [5] we show 
the results for four regions. The diagrams reveal that both the 
RGB and RC can be clearly discerned after the subtraction, con- 
firming that these features, and particularly the RC, are reliable 
detections. 
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Fig. 7. Top: CMDs of fields NW4, NW5, and NW6 in radial bins of width 10'. The radial angular distance from the M33 centre is 
displayed in the top-left corner of each plot. The isochrone that best fits the RGB slope and colour is also shown. The corresponding 
metallicity and age are indicated in the top-left corner of each panel. Bottom: CMDs in radial bins at angular distances greater than 
80' including stars in field NW7. The stellar populations in these regions appear to be older and more metal -poor with an average 
metallicity [M/H]= -1.3. 



4. The distance to M33 from the tip of the red giant 
branch 

Estimates of the distance to M33 obtained with different meth- 
ods vary considerably in the literature, ranging between 800 and 



940 kpc (|Lee et al.||2QQ2] |Ciardullo et aLl[2004| |McConnachie 
|et al.||2004[ |Tiede et al.||2004| |Sarajedini et al.|2006) . Here we 
derive the distance to M33 using the tip of the red giant branch 
(TRGB) as a standard candle ( jLee et al.|1993| l. The condition for 
the correct application of this technique is that the observed RGB 
should be well populated, with more than ~ 100 stars within 



1 mag of the tip (Madore & Freedman 1995 Bellazzini et al. 
|2001| l. Thus we built the / luminosity function of the RGB in 
the CMD of field SE3b, which is the one that most closely ful- 
fils this criterion, for all the stars redder than (V - 7)o > 0.7. 
We then applied an edge detection Sobel filter to the luminosity 
function to locate the position of the tip. The peak of the filter 
response is taken as the best estimate of the TRGB magnitude, 
while the half width at half maximum of the same peak (± 0.08 
mag) and the photometric error at this magnitude (± 0.03 mag) 
are taken as the associated uncertainties. An additional source 
of error is given by the uncertainty in the zero-point coefficients 



obtained from the photometric calibration, which in the I band 
is 0.05 mag. The tip is found at 7o = 20.73 ±0.1, with all errors 
combined in quadrature. Assumin g an absolute magn itude of the 
TRGB of Mj RCB = -4.02 + 0.05 ( jTiede et al.|2004}> , this corre- 
sponds to a distance modulus (m — M)Q = 24.75+0.11, implying 
a distance to M33 of 891 + 45 kpc. However, field SE3b contains 
stars in the farthest side of the optical disc of M33. As one can 
see from Fig. 6, the majority of the RGB stars used to determine 
the tip reside between 20 and 30 arcmin, corresponding to a pro- 
jected distance of 6 kpc or a deprojected radius of 10 kpc (being 
the field close to the minor axis). Thus the stars are about 8 kpc 
more distant than the center of M33, hence the implied distance 
to M33 is 883 ± 45 kpc. The uncertanties in the deprojection 
give an additional source of error of 2 kpc, which once added 
in quadrature, do not significantly increase the uncertainty in the 
distance. 

This distance modulus agrees with the results of ground- 
based observations of |Tiede etHl] ( |2004| ), (m - M) = 24.69 ± 
0.07, with HST observations of M3 3 disc stars, (m- M) () = 24.80 
+ 0.04 (random)+|J [^(systematic), dKim et al. 2002 1, and, within 
the errors, is compatible with the distance value of 840 + 40 kpc 
(Freedman et al.|2001| l that we assumed in Grossi et al. (2008). 
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Fig. 8. Mean colour of the RC as a function of the angular distance from the centre of M33 in the southern {top-left) and northern 
(top-right) fields. The size of the dots decreases with the distance to the M33 centre. A clear colour gradient is found in both regions. 
Mean theoretical RC colour and magnitude for different metal abundances (from [M/H]=-0.4 to [M/H]=-1.7) and ages are displayed 
in the bottom panels. The numbers along the model indicate different ages in Gyr units. The size of the filled dots in the lower panels 
is the same as in the upper panel to identify the different radial bins. 



5. Constraints on the age and metallicity of the 
stellar populations 

Evolved RGB and RC stars appear in a stellar population af- 
ter 0.9 -1.5 Gyr: stars first enter the RGB phase, where hydro- 
gen is burnt in a shell around the helium core, then move to the 
RC when the core-helium-burning stage begins. Even though the 
RGB and RC are affected by the age-metallicity degeneracy, we 
can use them together to constrain the mean metallicity and age 
of the stellar populations. As mentioned in Sect. 3.1, any varia- 
tions in the RGB colour and slope are more sensitive to metallic- 
ity; the magnitude of the RC has a stronger dependence on age 
( |Durrell et al.|2010) , while the colour of the RC for ages older 
than 2 Gyr is more sensitive to metallicity variations (Girardi & 
Salaris 2001). Bearing this in mind, we can use these features to 
perform a qualitative inspection of the CMDs at different radii to 
estimate the mean properties of the stellar populations as we pro- 
ceed to larger galactocentric distances. Thus, we built CMDs in 
radial bins of 10 arcmin width and compared the colour, the mag- 
nitude, and the shape of these features to synthetic isochrones 
and models to investigate possible radial gradients in age and 
metallicity. 

5.1. Age and metallicity gradients from the red giant branch 

Figures [6] and [7] show the spatial variation in the CMDs for the 
southern and northern fields, respectively. In each diagram, we 
compared the RGB feature within the magnitude range 20.7 mag 
< /<) < 24 mag to a set of stellar isochrones with metallic- 
ity between [M/H]=-1.7 and [M/H]=-0.4 ( |Marigo et ai:p 008). 
For each set of isochrones, we calculated the colour difference 
6(V - 1) between the observed RGB and the set of fiducial stellar 
tracks. The track that minimises the rms of the distribution of 
6(V - 1) is selected as an indicator of the mean age and metallic- 
ity of the stellar population in a radial bin. 



The average metallicity in the SE1-SE3 fields (Fig. [6]i 
changes from [M/H]=-0.4 at the edge of the optical disc (R < 
40', first two panels), to [M/H]=-0.7 at angular distances R > 
40', and the same abundance is obtained out to the outermost bin 
where a weak RGB is detected (70' < R < 80'). The correspond- 
ing mean age of the populations at large radii ( R > 40') does not 
vary considerably with the distance, being around log(age)=9.3- 
9.35. A blue plume of main-sequence stars is detected out to the 
third radial bin 40' < R < 50'. 

The properties of the stellar populations in the NW fields 
(Fig. [7]i appear similar to the SE ones when one considers the 
same range of angular radii (35' < R < 70'): the average metal- 
licity is found at [M/H]=-0.7 and the stellar ages in these ra- 
dial bins vary between 2 and 2.5 Gyr (9.3 and 9.4 in logarith- 
mic units). Young main-sequence stars are also detected in these 
fields out to a radial distance of 60'. The outer radial bins instead, 
at R > 70', show an older and more metal -poor RGB stellar pop- 
ulation, evidence that becomes more compelling when stars from 
NW7 are included (Fig. [7] Zower panel). The lower metallicity of 
the RGB stars at these radii is suggested by the overall bluer 
colour than the stars in the bins closer to the M33 centre. While 
the RGB and RC features become weaker from NW4 to NW6 
(upper panel of Fig. [TJ, there is a clear increase in the number 
of stars in NW7 out to the last radial bin (110' <R < 120'). The 
isochrones that most closely reproduce the slope and colour of 
the RGB have [M/H]=-l in the NW6 field (Fig.|7] upper panel), 
whereas when stars in NW7 are included (lower panel), the av- 
erage population at R > 90' appears to be even more metal-poor 
([M/H]= -1.3) and older (between 5 and 7 Gyr). 

5.2. Age and metallicity gradients from the red clump 

We derived the mean photometric properties of the RC as a func- 
tion of radius following Girardi & Salaris ( 2001). We determined 
the histogram of the / magnitude and colour distribution of stars 
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in regions of the CMDs including the locus of the RC, located 
within 24 mag < 7 < 24.7 mag and 0.7 mag < (V - 7) < 1.2. 
We then fitted a Gaussian plus a quadratic background function 
to both histograms, taking the central values of the Gaussians as 
the mean magnitude and colour of the clump. Given the diffi- 
culty in identifying a clear clumping of stars in the RC region 
in the third and fourth bin of the SE fields, or the low number 
of RC stars in the radial range between 70' and 90' in the NW 
ones, the results of our analysis at these radii should be treated 
with caution. To maximize the number of stars in the RC area, 
the range of angular distances between 70' and 90' was con- 
sidered as one single bin, and the outermost regions (R > 90') 
include detections from both NW6 and NW7. 

The results of our analysis are shown in Fig. [8] In the upper 
panels, we plot the RC colour versus radius for the SE (top-left), 
and NW (top-right) fields. The size of the dots decreases with 
the distance to the M33 centre to help us to identify the radial 
bins in the lower panel plots. In both regions, there is a clear 
change in (V - 7)o with radius. In the lower panels, the mean 
RC colour and magnitude are compared to the models of Girardi 
& Salaris (2001) to assess more reliably the range of age and 



[M/H]=-1.3 



abundances of the populations. The models show that variation 
in the metal abundances is the main driver for the colour change: 
the metallicity ranges from [M/H]=-0.4, [M/H]=-0.7 in the in- 
ner bins, down to [M/H]=-1.3/-1.7 in the outermost ones. This 
agrees with what we found from the RGB analysis. However, 
the error in the photometry (and in the distance modulus for Mj) 
displayed in the figure make the determination of the age subject 
to some uncertainties. In the SE fields, the redder bins corre- 
spond to stars with ages around 2 Gyr and [M/H]=-0.4, while at 
40' < R < 80' the colours and mean magnitude of the RC are 
compatible with an age between 0.8 Gyr and 1.5 Gyr, slightly 
lower than what we found for the RGB stars, and a metal abun- 
dance of [M/H]=-0.7. In the NW fields, the colour and magni- 
tude of the RC in the first bin may correspond to a population 
with either a mean age of ~ 2 Gyr and [M/H]=-0.7, or an age be- 
tween 1 Gyr and 2 Gyr with a higher abundance ([M/H]=-0.4). 
At 50'< R < 70', a metallicity of [M/H]=-0.7 would imply an 
age between 1 and 2 Gyr, although even older ages would be 
compatible with the photometric errors. 

At larger radii (R > 70') the NW data indicate that the prop- 
erties of the stellar population clearly differ. The colour becomes 
bluer than (V - 7)o = 0.85, corresponding to ages older than 5 
Gyr and an abundance of [M/H]=-1.3 (70'< R < 90') or lower in 
the most distant bins at R > 90' ([M/H]=-1.7). In the following 
section, we try to more tightly constrain the range of ages and 
metallicities of stars in the outermost field by comparing the ob- 
served CMDs to models with a range of star formation histories 
(SFHs). 

6. Star formation history of field NW7 

A more accurate constraint of the properties of the stellar pop- 
ulations in the northernmost region may help us to untangle the 
origin of the stellar feature around M33. Therefore, we applied 
the method of synthetic CMD fitting to the diagram of NW7 to 
derive the best-fit SFH of the stellar populations in this region. 

We used the StarFISH package ( [Harris & Zaritsky|[200T) , 
which determines the best-fit CMD according to a maximum 
likelihood statistic. To build the model CMDs, we selected four 
sets of isochrones with metallicities [M/H]=-1.7, -1.3, -1.0, -0.7, 
and ages between 1 and 10 Gyr (with a time binning step of 
log(age) = 0.2), as derived from the analysis of the RGB and 
RC features in Sect. 5. The code populates the set of stellar 
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Fig. 9. Star formation rate as a function of time of field NW7 ob- 
tained from the CMD fitting analysis (top), and the resulting syn- 
thetic CMD (bottom-left) compared to the observed one (bottom- 
right). A 6.3 Gyr (log(age) = 9.8) isochrone with [M/H]= -1.3 is 
overlaid on the RGB of both diagrams. 



isochrones (Marigo et al. 2008) assuming a given initial mass 



function (Salpeter 



1955), shifting the magnitudes according to 



the distance modulus, and applying photometric scatters derived 
from the artificial star tests (see Sect. 2.2 for details). Each syn- 
thetic CMD with a given range of ages and a fixed metallicity 
can be linearly combined to construct a model diagram for any 
arbitrary SFH. The code then determines the best-fit SFH per- 
forming^ 2 minimization of the differences in the number of stars 
between the model and observed CMD. 

We tried fitting several combinations of the four sets of 
isochrones and found the best-fit solution for a metallicity 
[M/H]=-1.3. The corresponding SFH and model CMD are dis- 
played in the upper and lower-left panel of Fig. [9] respectively. 
The error bars reflect the lcr confidence interval on the star for- 
mation rate averaged over an age bin. According to the recovered 
history of star formation, the stellar population in NW7 is metal- 
poor and older than 4 Gyr ago, and most of the stars formed 
around 6*}j Gyr ago. This agrees with our qualitative analy- 
sis in Sect. 5. The derived star formation rate corresponds to a 
total stellar mass of ~ 2 x 10 6 Ms . A comparison to the ob- 
served data (lower-right panel) shows a reasonably good match. 
An isochrone at 6.3 Gyr (log(age) = 9.8) is overlaid on both 
diagrams to guide the reader's eye. However, the match is not 
perfect since the colour of the simulated RGB appears roughly 
0.1 mag redder than the observed one at Iq > 23.5, and ~ 0.05 
mag bluer in the upper part of the RGB. Adding higher metal- 
licity isochrones at [M/H] = -1 to the model diagram in order to 
match the upper part of the RGB did not help us to improve the 
fit but led to an overall too red CMD. 
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Fig. 10. Spatial distribution of RGB and RC stars with 20.7 mag 
< I < 24.6 mag in the southern fields. The resolution of the map 
is l'x 1', smoothed with a gaussian kernel of FWHM five pixels 
(shown at the bottom-left comer of the figure). The colour scale 
indicates the number of stars per pixel. The stellar distribution is 
compared to the Hi complexes detected in this region. Contours 
are drawn at 3 to 10.5 XlO 18 cirT 2 in steps of 1.5 xlO 18 ctrT 2 . 

7. Stellar spatial distribution and comparison to the 
Hi features 



Figures 10 and 11 show the spatial distribution of the RGB and 
RC stars in the southern and northern fields, respectively. We 
selected stars with / magnitudes in the range 20.7 mag < I < 
24.6 mag, divided the fields into square bins with a size of l'x 1' 
, and determined the number of stars in each bin. The final image 
of the spatial distribution was obtained after convolving the 2-D 
histogram with a Gaussian filter of FWHM five times the bin 
size (5')- The density contours of the main Hi features detected 
in these regions are overlaid for comparison. 

The figures show several structures overlaid on a more dif- 
fuse distribution of stars. The main features in Fig.[T0|are: 

- The main overdensity at the top-right corner of the figure 
corresponds to the edge of the optical disc. A filament-like 
overdensity is connected to the stellar disc and extends out 
to the upper edge of the Hi complex (field SE3). 

- A fainter concentration of stars in the central part of the fig- 
ure (corresponding to field SE2) at a declination of 30°, pos- 
sibly related to the peak density of the main Hi complex. 

- An arc-like feature is visible to the southwest showing an 
enhancement in the distribution of stars in this region. This 
overdensity is not, however, spatially coincident with the 
southern extension of the Hi complex. It is possible that 
this increase in the number density of stars might be related 



instead to the stellar counterpart of the southern Hi warp 
HMcConnachie et al.|2009] >. 

In the northern fields, an extended stellar structure is detected 
in NW4 and NW5, along the Hi warp. Most of the stars appear 
to be located within the Hi disc, beyond which the stellar density 
rapidly decreases. Isolated and fainter clumps can be identified 
in field NW6. In particular, two overdensities match the posi- 
tion of one of the northern Hi clouds (whose column density 
contours are overlaid in the figure), although the stellar overden- 
sities appear to be offset from the peak of the Hi column density. 
Figure [TT] mainly stresses the difference in the stellar distribu- 
tion between the fields, along the Hi warp and the northernmost 
region. While the stellar density decreases from NW4 to NW6, 
field NW7 shows a higher concentration of stars, especially to 
the northern side of the field where few overdensities are visi- 
ble. The distribution of stars stretches to the north-west towards 
M3 1 forming a low surface brightness structure that extends to 
a maximum projected radius of ~ 3 degrees as one can see from 
the map of McConnachie et al. ( 2010 1. 



8. Discussion 

8.1. Hi clouds and stellar counterparts 

The search for optical counterparts to Hi clouds has often turned 
out to be a difficult task with very few successful detections. 
An optical study of 264 high velocity clouds down to a surface 
brightness limit of /xy ~ 26 mag arcsec~ 2 resulted in no detec- 
tions (Si mon & Blitz||2002| l. Even for one of the most massive 
high velocity clouds in the Milky Way, complex H (Hul sbosch| 
|1971| [Diet er 1971), infrared, and millimeter wavelength obser- 
vations did not reveal any evidence of either current star forma- 
tion or an evolved population ( |Simon et al.|2006] >. However, Hi 
clouds are often found close to stellar streams or disturbances 
in the stellar disc, as seems to be the case in M33. For exam- 
ple, the Hi clouds found near the southeastern edge of the disc 
of M31 ( |Thilker et al.1|2004l |Westmeier etaLl|20o31 |2008] l are 
spatially and kinematically correlated to the giant stellar stream 
(Ibata et al. 2001 ), implying that those clouds could be the re- 
sult of either tidal or ram-pressure stripping of a former satellite 
galaxy. 

The main purpose of our study has been to search for an op- 
tical counterpart to the most extended Hi features that we found 
in a 21 -cm survey around M33. The cloud properties seems to 
be indicative of a process of accretion or tidal disturbance in the 
outer regions of M33. This is confirmed by the features in the 
stellar population discussed in this paper and previous works, 
and also by the discovery of an extended cluster to the southeast 
of the galaxy (Stonkute et al. 2008). However we do not find 
convincing evidence of a young stellar population (~ 100 Myr) 
related to these features. This means that there has not been in 
situ star formation in these clouds, as indicated by their low col- 
umn densities, whose peak values in Hi are only around 10 19 
cm -2 (Grossi et al. 2008). We find a sparse and more evolved 
stellar population in the proximity of the Hi complexes in SE1- 
SE2, and some isolated clumps of stars are visible throughout the 
outer field NW6, but there is no clear evidence of a correlation 
between the stellar structures and the gas clouds. 

On the other hand, the stellar structures detected in fields 
NW4 and NW5 (Fig. [IT) follow the distribution of the Hi warp, 
suggesting that these stars represent the optical counterpart to 
the perturbed Hi disc, and that the stellar component - as in the 
case of the gaseous one - has a warped distribution. 
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Fig. 11. Spatial distribution of RGB and RC stars with 20.7 mag < I < 24.6 mag in the northern fields. The resolution of the map is 
l'x 1', smoothed with a Gaussian kernel of 5 pixels FWHM (shown at the top-right corner of the figure). The colour scale indicates 
the number of stars per pixel. The density contours of the Hi clouds detected in this region are overlaid as a comparison. Contours 
are at 3.5, 6, 8.5, 11, 13.5 xlO 18 cirT 2 . 



8.2. Stellar surface density at large radii: stellar halo or 
extended disc? 

The fields we observed with Subaru/Suprime-Cam probe dif- 
ferent regions beyond the optical disc of M33 (R25 = 35'. 5). 
Combining the information gathered from the analysis of the 
stellar populations, we found a trend of decreasing metallicity 
with radius. Determining the radial profile of the stellar surface 
density can help us to discriminate the nature of the different 
stellar components that are related to these metallicity variations 
at large radii (extended disc or a stellar halo). Thus, to derive the 
surface density we divided the surveyed area into radial bins us- 
ing two different geometries: a) circular rings to investigate the 
possibility that most of the stars we detect trace a smooth stel- 
lar halo; b) elliptical annuli with position angle and inclination 
given by Corb elli & Salucci"| ( |2000| l, if most of the stars lie in an 
extended disc that follows the Hi distribution. We discuss these 
possibilities in the following subsections. 

We selected RGB stars with Iq between 20.7 mag and 24 
mag, and colours (V - 7)o > 0.7 mag, limiting our analysis to 
RGB stars to ensure a higher completeness level on the photom- 
etry. 

8.2.1. Stellar halo density profile 

We derived the number counts in radial circular bins of width 
5' for R < 40' and width 10' for R > 40'. A different spac- 
ing was chosen to more accurately trace the profile at the edge 



of the M33 stellar disc. The number counts in each bin were 
then divided by the area of the corresponding region. The level 
of contamination was estimated from the control field selecting 
the number of detections within the same magnitude and colour 
range of the RGB stars (see Fig. |4j>. The surface density of con- 
taminants was found to be ~ 0. 1 arcmin -2 , and subtracted from 
the number counts in each radial bin. 

Figure [12] {left panel) shows the RGB surface density versus 
the projected radius. The open squares and the filled dots indi- 
cate the contamination-subtracted surface density in SE1-SE3 
and NW4-NW6, respectively. Star symbols represent the sur- 
face density when NW7 is also taken into account. Error bars 
are given by the square root of involved star counts scaled by the 
area of each region. 

Our data set allows us to study the decline of the M33 den- 
sity profile out to a projected radius of 90' and 120' in the SE 
and NW fields, respectively. A change in the shape of the profile 
occurs in the SE fields (open squares) at radii greater than 40', 
where the optical disc ends; beyond 50', the profile shows a sim- 
ilar trend in both the northern and southern regions. The density 
enhancement at radii greater than 90' (star symbols) is due to the 
population in field NW7 belonging to the large stellar feature to 
the north of M33. The density here is about two times higher 
than that found in NW6 only. The overdensity extends further 
north of NW7 for about one degree along the direction to M31 
(McConnachie et al. 2010). 

A power-law profile, R", representing the spatial distribution 
of a spheroidal stellar halo, was then fitted to the data at R > 40' , 
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Fig. 12. Stellar surface density profiles in the observed fields versus projected (left) and deprojected (right) radius. The open squares 
and filled dots show the stellar surface density in the SE1-SE3, and NW4-NW6 fields, respectively. Star symbols indicate the density 
at large radii when field NW7 is also included. The M33 stellar disc extends out to 40', while at larger radii there is a clear change in 
the surface density. A power-law profile (R") is fitted to the data in the left panel to investigate the presence of a stellar halo. We did 
not attempt to fit the star points because of the anomalous overdensity in field NW7. An extended disc component (e~ R ^ h ) is fitted 
to the data in the right panel. The best-fit parameters for each model are displayed in the panels. 



leaving out the overdensity in NW7 (star symbols). We find best- 
fit power-law indices n = -2.13 ± 0.24 and n = -2.24 + 0.33 
in the NW and SE regions, with ax 2 of 1.4 and 1.8 , respec- 
tively. The two resulting power-law indices are compatible with 
a unique distribution, even though the northern side has a some- 
what higher density. Our result is comparable to that found in 
the MW and M31 stellar halos: the surface density profile of 
halo stars in the Galaxy falls off with galactic radius following a 
power- law index varying between -1.5 and -2.6 (Chiba & Beers 
2000, |Juric et al.||2U08l |Bell et al.||2008l |Juric et al ||2008| l. In 
M31, Ibata et al. (2007) obtained a surface density profile with a 
power-law index of n = -1.91 ±0.12, while Tan aka et al.| ( |2010] ) 
found a steeper power law with n = -2.17 ± 0.15. 

8.2.2. An extended outer disc 

We investigate the presence of an extended stellar disc under the 
assumption that it follows the same orientation as the warped Hi 
disc. If an extended stellar disc is in place, the stellar surface 
density can be inferred by deprojecting the stars onto the plane 
of the galaxy. For a warped disc, a tilted ring model is needed. 
Following [Corbelli & Salucci| ( |2000| l, we assumed that the outer 
disc position angle varies between 21° and -11°, while its incli- 
nation stays roughly constant at about 50°. Hence, we computed 
the location of the elliptical contours in the plane of the sky cor- 
responding to fixed galactocentric radii, i.e. the radii of the tilted 
rings. We used radial spacing of 5' for the innermost ellipses and 
10' for ellipses with major axis larger than 50'. We then assigned 
each star to a ring and compute the stellar surface density in the 
ring areas covered by our survey. 

The deprojected surface density is displayed in the right 
panel of Fig. [12] The figure now shows an excess density in the 
SE fields between 40' and 60' , and a very good match between 
the density values of the northern and southern regions in the 
following two bins (60' < Rj ep < 80'). An exponential profile 
was fitted to the density distribution for Rdep > 40' . The best-fit 
scale-length is h =25'.4 ± 0'.3, with^ 2 = 0.7. This scale-length 
(~ 7 kpc) is much larger than the stellar scale-length of the in- 



ner disc (~ 1.5 kpc). The stellar drop-off in M33 is also slightly 
shallower than the Hi decline in the same area. 

Given the better fit obtained after deprojecting the stellar dis- 
tribution, we conclude that the observed surface density in our 
fields seems to be more consistent with an extended disc compo- 
nent than a smooth stellar halo. However, we emphasize that our 
analysis is based on fields selected along the Hi warp, and our 
results do not exclude that the contribution of a faint spheroidal 
halo could be investigated more effectively in different regions 
around M33. Fields from the PAndAS survey at similar angu- 
lar distance as NW6 and NW7 along directions not polluted by 
either the disc or the large stellar feature around the galaxy, con- 
tain a small RGB population and are presumed to represent an 
underlying halo population in this galaxy. 

8.3. The large stellar feature around M33 

The origin of the large stellar feature around M33 is discussed 
in M cConnachie et ak] ( |2009| |2010| l, and we refer the reader to 
those papers for more details. Here we summarise the main con- 
clusions of these authors in light of our additional results. Their 
favoured scenario assumes that the structure is the result of a 
tidal disturbance caused by the motion of M33 around M31. The 
northwest - southeast symmetry of the "S-shaped" structure, and 
the rough alignment with the Hi warp, are the main pieces of 
evidence supporting this interpretation. Simulations of the M31- 
M33 system have showed that a close encounter at a pericentre 
distance of about 40 kpc would excite tidal tails in M33 with- 
out severely distorting or disrupting the disc. These simulations 
also set the epoch of the encounter to be between 2 and 3 Gyr 
ago. However, within this scenario it would be difficult to rec- 
oncile the metallicity difference between the structure ([M/H]=- 
1.3) and the disc ([M/H] =-0.7,-0.4) stellar populations, a differ- 
ence that we also detect within our data set. 

On the other hand, if the structure originated from the disrup- 
tion of a dwarf galaxy such a difference in the metal abundance 
would not be an issue, and this scenario is also considered as a 
plausible alternative. The lack of a visible progenitor would im- 
ply that the dwarf galaxy has been entirely destroyed prior to the 
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Fig. 13. Hi surface density versus deprojected radius. The verti- 
cal line shows the maximum distance at which young blue stars 
were detected. Young stars are found out to a radius that is com- 
parable to the distance beyond which the Hi surface density de- 
clines abruptly. 



current epoch. From the SFH that we derived in Sect. 6, the to- 
tal st ellar mass in field NW7 w ould amount to ~ 2 x 10 6 M Q , 
and |McConnachie et al. ( 20T0] > estimate that the total luminosity 
of the whole structure is My = -12.7, comparable to that of a 
bright LG dSph. Our data cannot provide constraints that enable 
us to rule out one of these possibilities, but it is interesting to 
note that the lower limit we infer to the age of the stars (log(age) 
= 9.5 or ~ 3.2 Gyr), is comparable to the epoch of the encounter 
between the two galaxies predicted by the simulations. 



8.4. Recent star formation beyond R25 

In Sect. 3, we pointed out that a population of young stars, with 
ages around or older than 200 Myr, was detected in SE3, NW4, 
and NW5. The majority of them were found within 50', with a 
few blue stars detected out to a galactocentric radii of 60'(~ 17 
kpc; see Fig. |6]and|7]i. Finding evidence of recent star formation 
activity is surprising at these radii, given the low column density 
of the gas. As a comparison, we show in Fig. 
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_ the M33 Hi 

surface density as a function of the galactocentric distance along 
the Hi warp. At 60', the gas surface density is around 0.6 M Q pc~ 2 
(or ~ 7.5 xlO 19 cm -2 ), roughly eight times lower than the density 
at R25, then starts to decline abruptly at Rd ep > 70'. 

The presence of such young populations at large galactocen- 
tric radii challenges the current picture of star formation given 
by the Schmidt-Kennicut law ( |Kennicutt|1998| l. However, stud- 
ies of nearby spirals based mainly on GALEX observations have 
shown that star formation can also occur beyond ^25 in a low- 
density environment (Thilker et al. 2005a, 2007 ). Stellar asso- 
ciations with an average age of ~ 150 Myr were found in the 
outer discs of a sample of five nearby spiral galaxies (Alberts 
et al. 2011). |B~igiel et aL] ( |2010[ ) studied star formation rates in 
the far outer discs of spirals where Hi is the dominant gaseous 
component. In contrast to the inner parts of the disc, the Hi den- 
sity and the star formation rate (derived by the far UV emission) 
showed an overall good correlation, implying that different pro- 
cesses may govern the formation of stars at large galactocentric 
distances. 



9. Summary 

We have presented Subaru/Suprime-Cam V and / observations 
of seven fields around the disc of M33. The fields cover regions 
of the sky where Hi complexes were detected by the ALFALFA 
survey, the northern Hi warp, and part of the extended feature 
that was discovered around the galaxy by the PAndAS survey. 
We have studied the properties of the different stellar populations 
using CMDs, estimated the mean metallicity and age in each 
region, and determined the radial distribution of stars. Here we 
summarise our main results: 

1. We confirmed the presence of an extended population of 
evolved stars (in the RGB and RC phase) with ages between 
1 and 10 Gyr out to an angular distance of 120' (~ 30 kpc). 

2. A sparse population of young stars with ages around 100- 
200 Myr is detected both to the NW and to the SE out to a 
galactocentric radius of ~15 kpc. 

3. We did not find a clear optical counterpart to the main Hi 
complexes detected around the optical disc, although stellar 
structures are detected in all target regions. This provides ev- 
idence that the stellar distribution at large radii is disturbed 
in a similar way to the gaseous one because of either tidal 
interaction or accretion. 

4. We used the shape and colour of the RGB as well as the mag- 
nitude and colour of the RC to derive the metallicity of the 
dominant stellar population in each field. The mean metallic- 
ity gradually decreases from [M/H]=-0.4 at the edge of the 
optical disc, down to < -1, out to a radius of ~ 30 kpc. 

5. The surface radial profile of the stellar distribution in our 
fields can be more closely fitted by an extended disc com- 
ponent with a scale-length of ~ 7 kpc, rather than a smooth 
stellar halo. This suggests that the stars as in the case of the 
gas have a warped distribution. However, this does not ex- 
clude the presence of an underlying stellar halo, which may 
be more accurately studied by targeting other regions around 
the galaxy. 

6. Finally we analysed the stellar populations in a region be- 
tween 1.5 and 2 degrees north of M33 probing part of the 
large stellar plume surrounding the galaxy. The RGB and RC 
colours show a more metal-poor and older population than in 
the other fields. A CMD fitting analysis shows that stars in 
this region formed between 3 and 10 Gyr, with a metallicity 
of [M/H]=-1.3. 
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